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Abstract 

In this article, we investigate the possibility of approximating the physical inner 
product of constrained quantum theories. In particular, we calculate the physical 
inner product of a simple cosmological model in two ways: Firstly, we compute it 
analytically via a trick, secondly, we use the complexifier coherent states to approx- 
imate the physical inner product defined by the master constraint of the system. 
We will find that the approximation is able to recover the analytic solution of the 
problem, which solidifies hopes that coherent states will help to approximate solu- 
tions of more complicated theories, like loop quantum gravity. 



1 Introduction 

Loop Quantum Gravity is an attempt to quantize General Relativity to obtain a 
quantum theory of gravity (see [19] and references therein). In order to do so, one 
first has to rewrite GR as a Hamiltonian theory, which leads to a constrained sys- 
tem. The question how to quantize a Hamiltonian theory with constraints exists 
for a long time [10, 6, 4]. The technical issues are, in fact, not too complicated. 
The constraint functions Cj, which are phase space functions in the classical the- 
ory, become operators Cj on an auxiliary, or kinematical Hilbert space TCkin- The 
classical issue of solving the constraints is, on the quantum side, transformed into 
the task of computing the common zero of the constraint operators Ciip = 0. Un- 
fortunately, since the point zero is generically in the continuous parts of the spectra 
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of the Ci, the space of solutions to the constraints will in general be no subspace 
of the Hilbert space. Rather, the set of solutions will contain distributions, which 
will not be normalizable in the inner product on TCkin- So, even when the space of 
solutions of Ciip = is known, it has no structure of a Hilbert space yet. In fact, 
from various examples it is known that the choice of a " physical" inner product is a 
quite nontrivial task, and the resulting theory will crucially depend on this choice. 
Generically, a lot of choices are possible which result in trivial theories, or nontrivial 
theories that in no way resemble the features of the classical theories one started 
with. 

In Loop Quantum Gravity this issue is even more complicated, since the alge- 
bra of constraints is tremendously difficult. Although it is possible to define the 
Hamiltonian constraints [13, 14] (which then still contain a lot of freedom), the 
diffeomorphism constraints cannot be implemented as operators on Hkin [19]. Fur- 
thermore, there are arguments that even if such operators could be constructed, 
they could not be self-adjoint, which makes spectral analysis quite difficult (see, e.g. 
remarks in [19, 12]). Nevertheless, a diffeomorphism invariant Hilbert space Ti-diff 
could be defined, which contained vectors invariant under the exponentiated diffeo- 
morphism operators [3]. But since the Hamiltonian constraints do not even weakly 
commute with the diffeomorphism constraints, there was no way of defining the 
Hamiltonian constraints on TCdiff, so the search for the solutions to the constraints 
could not continue. 

In 2003, a solution to this set of problems was suggested [12]. It is mainly mo- 
tivated by the observation that a classical constrained system with an arbitrary 
number of constraints Cj can be replaced by an equivalent system with only one 
constraint, the so-called master constraint M. Solving for all the Ci simultaneously 
or solving for M then leads to identical solutions. In the corresponding quantum 
theories the constraint operators Ci can be replaced by the master constraint oper- 
ator M. Unfortunately, while on the classical side these two sets of constraints lead 
to equivalent classical theories, this equivalence is not that obvious in the quantum 
theory. In generical situations, there are a lot more solutions to Mtp = than to 
Citp = for all i. On the other hand, in generic models, even not all of the solutions 
to the constraint, whenever one was able to compute them, had non- vanishing phys- 
ical norm, electromagnetism being the first example. In many tested models, many 
of the formal solutions of M = also turned out to have zero norm in the physi- 
cal inner product [5]. This reinforced the hope that the master constraint could be 
used to define the physical inner product and lead to sensible results: There was the 
possibility that the additional solutions were all spurious, i.e. had all zero physical 
norm. 

Still, even with the master constraint, the task of computing the physical inner 
product is, for the case of Loop Quantum Gravity, technically very hard [16]. The 
operator M is quite complicated and hence a spectral analysis is presumably not 
possible analytically. On the other hand, this is nothing to be severely concerned 
about, since the knowledge of the physical inner product would be equivalent to a 
complete solution of the theory, i.e. a quantum theory which, in its semiclassical 
sector, is believed to contain classical general relativity. But even of classical GR 
one does not know all solutions (in fact, one does not even know how many solutions 
one doesn't know). So solving the complete theory is not something we can expect 
at all. On the other hand, this is also not something that is necessary in order to do 
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physics. The situation is similar to the case of condensed matter physics, where the 
Hamiltonian operator for a crystal with roughly 10 23 particles can be written down, 
but not one solution is known. And this is not necessary in order to do physics, since 
no exact solutions are needed to make predictions about conductivity, magnetization 
or the shear module of a crystal. 

So what one needs in order to extract physics from the equations of LQG, are 
approximation schemes for the physical inner product defined by M. In this paper, 
we will investigate a simple cosmological model with only one constraint C and 
try to approximate the physical inner product defined by the corresponding master 
constraint M with the help of complexifier coherent states. These states have been 
introduced by Hall [7, 8, 9] and investigated in [17, 20, 21]. They have properties 
that make them useful tools for semiclassical analysis, in particular, they will allow 
for an approximation for the physical inner product of our toy model. 

There are, with some tricks, other possibilities to calculate the physical inner 
product in this model analytically [11], without the need to employ the master 
constraint M, and we will see that the approximation provided by the coherent 
states in fact works. Furthermore, the solutions to the master constraint M that 
are no solutions to the constraint C turn out not to be in the physical Hilbert space, 
as one hopes. 

2 Preliminaries 

2.1 The master constraint 

Consider a phase space M. and a set of first class constraints Q. This means that 
the Poisson brackets between two constraint functions {Q, Cj} vanish at the phase 
space points where all the constraints Q vanish. Then the space of all points m <E M 
with Ci(m) = for all i is called reduced phase space. The Hamiltonian vector fields 
of the constraints Xq generate flows that leave this reduced phase space invariant, 
and the orbit space is called the physical phase space. All phase space functions 
weakly commuting with the constraints, called weak Dirac observables, can then be 
made into functions on the physical phase space, and will serve as observables. 

Instead of solving all the (possibly infinitely many) equations Ci{m) = 0, one 
can also define the so-called master constraint [12, 16] 

M := ^CiKijd. (2.1) 

i 

Here Kij is a symmetric, positive definite matrix in the case of i being a discrete 
set. Otherwise, has to be a positive definite operator kernel and the summation 
over i turns into an integration. It is straightforward to see that 

M(m) =0 4* d(m)=0 for all i. (2.2) 

Furthermore, for any phase space function / weakly commuting with the constraints: 

{ {M, /}, /} « 0, & {d, /} « for all i. (2.3) 

So M enables us to derive the complete set of observables on the physical pase 
space. This means that the physical phase space itself can be constructed from the 
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knowledge of M, so one does not lose any information if one goes over to M from 
the Cj. The final classical systems defined by both are in fact equivalent. This is 
in fact independent of the actual choice of Kij, so there are a priori many possible 
master constraints. One can choose the one that is most useful, makes the sum 
(2.2) converge and is the most convenient to compute. This freedom is quite useful 
in the quantized theories [5] . 



2.2 Quantizing constrained systems 

Consider a Hilbert space TCkin an d a se t of constraint operators Cj. Let furthermore 
S>kin C Ti-kin be a dense subspace which is invariant under the action of the Cj and 
their adjoints c\. A solution to the constraints is a linear form 



-»• C (2.4) 
with 

I o C\ = for all i (2.5) 
The space of all such linear forms @* hys is a subspace of the algebraic dual of S^kin'- 

%hys C % m . (2.6) 

In general, ^^ ys H TCkin = {0}, so the solutions to the constraints will not be 
normalizable in the inner product on TCkin- Furthermore, not only can S>* hys not 
inherit the inner product from TCkin-, experience also shows that not all elements in 
S>* hys are acceptable as solutions for generical models. Generally, S>* hys is "too big" 
to admit a positive definite inner product. 

The task is now to find a Hilbert space TCpfi ys that is a subspace of ^p^yg- 

In 

general, there will be a lot of them, and the issue of finding the "right" one is a 
nontrivial question. One necessary condition for this Hilbert space would for in- 
stance be that it supports the quantizations of observables weakly commuting with 
the constraints Cj as self-adjoint operators. 

The idea is now to replace the set of constraint operators Cj by the quantization of 
(2.1), the master constraint operator 

M := ^CiKijCl (2.7) 

i 

Here the freedom to choose has in fact to be exploited in order to arrive at a 
well-defined operator. 

One can immediately see that 3>kin is left invariant by M, and for any linear 
form I e ®* kin : 

I o Cj = for alU I o M = 0. (2.8) 

If is only contained in the discrete (point) spectrum of all the Cj, then (2.8) is 
an equivalence relation. But if is in the continuous part of the spectrum, then 
the converse of then (2.8) need not be true. So the space of solutions of M is 
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potentially larger than the space of solutions to the constraints Q. On the other 
hand, M provides a natural way to define the physical Hilbert space. In all models 
checked so far, this Hilbert space turned out to contain no solutions of M that were 
no solutions to the C{. Although there are no hard proofs available at the moment, 
there are high hopes that this is the case generally. 

2.3 The physical Hilbert space and rigging maps 

Since M is formally symmetric and non-negative by construction, one can always 
extend it to a self-adjoint operator. Of course, this has to be checked from case 
to case, and we assume it to be true here. Then M is also positive and hence its 
spectrum is a subset of [0, oo). We furthermore assume that is in fact in the 
spectrum of M. Then, for each A € spec(M) there is a Hilbert space H\, equipped 
with inner product ( • | • )\. Furthermore there is a measure /x on spec(M) such that 
each vector of TLkm can be decomposed into 

V> = / _ dji(A) Y>(A) (2.9) 

J spec(M) 

and the inner product in Ttkin via 

<V#> = / . ^(A)<V(A)|0(A)> A . (2.10) 

J spec(M) 

These facts can be abbreviated by 




(2.11) 



which is called the direct integral decomposition of TCkin- We neglect all discussions 
about the mathematical subtleties of this construction and just state the fact that 
the Hilbert space 7i\ for A = shall serve as physical Hilbert space. 

Hphys '■= Wo, (-\-)phys '■= ( ' I ' )o- (2-12) 

This procedure can in principle be applied to every system with constraints. In 
particular, in the case of being in the discrete part of the spectrum only, this is 
just the definition of zero eigenspace of M. In the case of a mixed spectrum of M, 
one has to decompose Hkin a "discrete" and "continuous" part and perform 
this analysis for each part separately: 

n kin = nl p m © ntn = [ e . d»pp(\) nf © r . cWc(A) Hf. 

JspeCpp(M) J spec ac (M) 

(2.13) 

Here "pp" and "ac" denote the pure-point- and the absolutely continuous part of 
the spectrum of M. In the end, the physical solutions from each separate part 
constitute the whole physical Hilbert space, in particular 

n phys = n p p © n% c . (2.14) 
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Note that Hq P C Hkin consists of the zero eigenvectors of M. For both cases "pp" 
and "ac" the Hilbert space Ho is a subspace of ^t^ys m ^he following sense: In 
the direct integral decomposition ip £ Hkin corresponds to a map A i— ► ip{\) with 
■0(A) £ H\. Denote the evaluation of this map at A = as tp(0). Then define for 
each ip £ Hkin the following linear functional: 

: <p (m\m)o- ( 2 - 15 ) 

So rj provides a "projection" of Hkin i n Ho (which is a proper projection for the 
pure-point case). On the other hand, since the decomposition of M<p is given by 
the map A i— > A</>(A), it is immediately clear that rj as defined in (2.15) maps ip into 
a solution of the constraint M, i.e. rj[ip] o M = 0. In the literature one also finds for 
rj the name rigging map. 

Although conceptually quite clear, the decomposition (2.11) is in general nontrivial, 
as the spectral analysis for an arbitrary self-adjoint operator can be highly compli- 
cated. This is also the fact for the master constraint operator of LQG. Although 
the operator has been constructed explicitly [16], its complete spectral analysis is 
currently out of reach. 



2.4 The operator 5(M) 

We now would like to find a way of approximating the solutions of M, i.e. the linear 
functional (2.15) for all ip € Hkin- For this we observe 

7?[v# = <V(0)|^(0)> 

= / . dji(X) <V(A)|<5(A)|<KA)) A (2.16) 
= (i(,\8(M)\<f>), 

where S(M) depends on dji(X). The "operator" S(M) is, of course, no operator but 
a quadratic form which will not exist for all tp, <p £ Hkin-, but only for those vectors, 
whose decomposition into A i-> ^(A) is regular at A = 0. 

Equation (2.16) now enables us to express the linear form r)[ip] without directly 
computing the direct integral decomposition. All we need is the construction of a 
(5-function for the measure fi on the spectrum of M. For this, the knowledge of the 
measure jx is, of course, sufficient, but not necessary. All one needs to know is the 
rate at which the measure "diverges" or "goes to zero" at A = 0. 

If is in the discrete part of the spectrum, the measure ji is a Dirac-measure 
there, and the delta-function is just a usual Kronecker-delta. In this case the integral 
(2.11) is just a sum over A = and hence <5(M) is a proper projection operator which 
maps Hkin onto its subspace Hq P . 

If is in the continuous part of the spectrum, the measure fi can behave very 
singularly at A = 0. This is supported by the following observation: Consider the 
didirect integral decomposition of a self-adjoint operator A, whose spectrum is the 
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real line, and whose spectral measure fi^ is the Lebesgue measure dX. Then the 
operator A 2 has [0, oo) as spectrum, and its measure fi^ 2 is then -^%dX, which 
diverges as A tends to zero. So by taking powers of operators, one changes the 
behavior of the spectral measure at A = 0. 

So, an educated guess for a delta-sequence are the functions 
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A .— ► r a e-^r. (2.17) 

If the parameter a £ M. is chosen correctly, so as to match the degree of divergence 
of the measure [i at A = 0, then this diverges to a multiple of the ^-distribution at 
A = 0, as t tends to 0. Determining a can be nontrivial. On the other hand, there is 
a different way of calculating the physical inner product without the need to know a 
explicitly: Choose a vector ij) G Hkin such that in its direct integral decomposition 
one has 

< ||^(0)|| < oo. (2.18) 

Then one has 

M 2 _ _m£ _ M 2 

e t t a e * t e * 
lim = lim = lim 

^° We-^W) ^° t-« <V|e-^|V> ^° /d/x(A)(V(A)| e -^|V(A))A 

lim t _> t~ a e 



lim^o s -« Jd/x(A)(V(A)|e--|V(A))A 
1 



(2.19) 



k(0)|| 



<5(M). 



Here, the equalities have to be understood in the sense of quadratic forms. This 
shows that, by employing a "reference vector" ip one can save oneself from the need 
to explicitly calculating a. The resulting physical inner product is then rescaled 
such that r][ip] has unit norm in the physical inner product. One now has shifted 
the problem of finding a to finding a reference vector that satisfies (2.18). But 
since the set of vectors satisfying (2.18) is dense in TLkin for A = being in the the 
continuous spectrum of M, one would think that finding such a vector is not too 
complicated. 



3 The model 

As a test of full LQG one can focus attention to symmetry truncations of LQG, 
in the hope to get technical insight in soluable situations which are still valid in 
the full theory. A particular class of truncations modelling cosmological situations, 
called Loop Quantum Cosmology (LQC), has received a lot of interest recently 
[2]. We will apply the approximation scheme we have in mind to the easiest LQC- 
model possible: that of a spatially homogenous and isotropic universe without any 
matter fields. For the case of a nonzero cosmological constant, the physical Hilbert 
space has been constructed in [11]. The case for vanishing cosmological constant is 
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different, since then the zero is not in the discrete part of the spectrum any longer. 
It is this case that we will investigate in this article. 
The metric the model describes is 

ds 2 = -dt 2 + a{tf {dx 2 + dy 2 + dz 2 ) 

The canonical variables that resemble the Ashtekar variables in full LQC are 

P = a\ c = A, p.2) 

which have 

{c,p} = 1. (3.3) 

In the case of Loop Quantum Cosmology, p is allowed to actually take values in R, 
so as to allow for a change in orientation. Note that this normalization differs from 
the literature by a factor of Stt^Gn /3 [2]. 
The constraint of the system is given by 

C = -24sgn(p) • ^pjc 2 + 3Av^. (3.4) 

The kinematical Hilbert space TCkin of LQC is the space of almost periodic functions 
in c, which has the functions 

\ v ) ■= c i — ► e*5 c , is£R (3.5) 

as a basis. So TCkin is n °t separable, since it has an uncountable basis, a feature that 
it shares with the Hilbert space of full Loop Quantum Gravity. The quantization of 
p acts on TCkin via differentiation by c: 

PW) = ^1"), (3-6) 

whereas c is not definable as an operator. Rather, one has to employ exp(zAc) for 
all A € R as multiplication with the corresponding function of c, via: 

e iX£ \u) = \\ + 2v), (3.7) 

which shows that the failing weak continuity of A ^ exp(iAc) is the reason why 
one cannot define c without the exponential. This amounts to a problem when 
attempting to define the quantization of (3.4), since it explicitly contains c. This 
problem occurs in a similar fashion in the full theory, and it is resolved in LQC by 
introducing a parameter no and approximating c 2 by /Xq 2 sm 2 (noc). 1 

If c <C l//Uo, the two functions are quite similar, thus, if is sufficiently small, 
the function will be a good substitute for a large part of phase space. The lat- 
ter function, however, can be quantized on TCkin-, which results in the following 
constraint operator: 

C = -24 ( ^ V sgn p ^ + 3A^ § , (3.8) 

1 In full LQG the parameter no is not necessary, because spatial diffcomorphism invariance "swallows 
it up" there. In LQC the spatial diffeomorphism group has been gauge-fixed so that one has to deal with 
fio as an artifact of the model 



(3.1) 
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where 



sgnp VIpI W) = o \ v ) ( 3 - 9 ) 

6fl 

is the quantization of sgn p \/jp[, that results from the attempt to keep the quanti- 
zation scheme as close as possible to the ones employed in Loop Quantum Gravity. 
The introduction of the scale induces a splitting: 

n km = H 5 (3.10) 

<Se[o,4no) 

where each Tig = span{|z^), v € S + 4^o^} is separable and left invariant by the 
constraint. Thus, one can calculate the rigging map for each sector 5 € [0,4/xo) 
separately. 



3.1 Solutions to the Constraint 

We attempt to calculate the physical inner product for the above cosmological model 
with A = 0. For this, we first use a trick that has been employed in [11] to compute 
the physical inner product for the case with A / of Riemannian cosmology. 
For this case the is in the pure point spectrum of the constraint, and hence the 
eigenvalue problem could be solved explicitly. It was reasoned that this stays true 
for the case of Lorentzian cosmology with nonzero cosmological constant. If A is 
set to zero, then Riemannian and Lorentzian cosmology become indistinguishable, 
and then is not any longer in the pure point spectrum of C. Still, the trick used 
in [11] will make it possible to compute the physical Hilbert space analytically. 

Since the trick is in principle not applicable in full Loop Quantum Gravity, we 
will also try to calculate the physical inner product without using it. Furthermore, 
since without the trick the constraint operator is no longer symmetric, we work 
with the master constraint operator M = C^C instead. In this setting, the master 
constraint operator will be too complicated to derive the physical inner product 
analytically, but the coherent states will enable us to approximate it well enough. 
This is precisely what one expects also in full LQG In particular, we will find the 
approximate solutions being close to the solutions previously obtained analytically, 
which will show that the coherent states are in principle useful to approximate phys- 
ical inner products. It will also demonstrate that, although the original constraint 
operator C has been replaced by M, which has potentially many more solutions, 
the physical Hilbert space computed from M consists of solutions of C only. 

3.2 Riemannian cosmology with and without cosmolog- 
ical constant 

We shortly repeat the results of [11]: There, the Hilbert space was the same that 
the one used in the Lorentzian case we are looking at, but the constraint looked 
differently. In particular, it was given classically by 

a 

CRiemann = 24sgnp-^| C 2 + 3AyW. (3.11) 
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The trick is to modify this constraint by multiplication with \p\ 2 , to arrive at the 
(up to the points p = 0) classically equivalent function 

CRiemann = 24 C 2 + 3Ap. (3.12) 

By introducing a length parameter [1q, as shown in section (3), this function can be 
quantized on Ttkin as 

CRiemann = 24 f^) 2 + 3Ap, (3.13) 

V w / 

from which the physical inner product was computed in [11]. The only solution to 
(3.12), which has nonzero physical norm, was calculated to be 

I = V J +-^ (t^-s) (A|, (3.14) 



*M0 2A^f. \2A/iX / 
A M 2 

where J n (x) is the n— th Bessel function of first kind. In particular, although there 
exist solutions to the constraint in every sector 5 £ [0,4^o)> only the solution in the 
sector 5 = 4/xo[l — mod ( A 3 3 , 1)] has nonzero physical norm. Here, mod (x, y) 

denotes the remainder of x/y. So by the constraint, one sector is allowed, and the 
resulting physical Hilbert space is one-dimensional, as assumed for a theory with 
zero degrees of freedom. 

If one lets A — > 0, neither does the solution (3.14) converge, nor does the sector 
5 in which it is defined. So, in order to investigate the case A = 0, we have to 
calculate the physical inner product explicitly, instead of taking the limit A — ► of 
the already known solutions. Classically, this amounts to 

C = 24c 2 . (3.15) 

Since the quantized constraint 

^ = 24 sinVoc 

is self-adjoint, we do not employ the master constraint, but calculate the rigging 
map (2.19) directly with C. We have to compute (2.19) 

c 2 

(n\ e~— \u) 



v[W)]\fj) = ]jm _c2 

~~* {v | e t | u ) 

where \vq) is a reference vector for correct normalization. 
Since C is bounded, we have 

?1=0 

2ra 



^ (-1)" / 24 V" , , . 4 „, „, v 
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With definition of e 4W)C one can see that 
(/x|sin 4 "(/x c)|i/) 



f v—\x / 4??< \ 

4»,...*,..x _ J M~V** ( 2n+ ^ ) for^ E Z (317) 

else 



So in what follows, assume to be integer. Then we have 
, , 62 , v , „a ^ ("l) n / 6 \ 2n / 4n \ 

l) n / 6 \ 2ri r(4n + l) 



, ^ (-l) n / 6 \ 



to n] U^/ r(2n + l + ^)r(2n + l-^) 

- - r _iyv v izil! r(4n + 1) 

- <> l > 2_. n [ T{2n + l + N)T(2n + l-N) ^ 8j 

n=0 

with the abbreviations 



^ v — n 6 
-/V = — , a - 



Now we use the identity 

n=0 

(see e.g. [1], p. 559), for / being as in (3.18). A coordinate transformation s = 
yields 



£ 

n=0 



("l) n a 2n r(4n + l) 



n! r(2n + l + A^)r(2n + l-iV) 



1 /*ioo 

= — / dsr(-s) 

2vri 7_ loo 



r(4, + 1) 



r(2s + i + iv)r(2s + i-iv) 



2 



8tt* y i00+1 r(i-| + Ar) rQ-|-iv) 

Using again (3.19), we get 



E 



(-l)" a 2n T(4n+1) 



n! r(2n+l + iV)r(2n + l- N) 



1 

4 a 2 g n! " " r(i-f + iV)r(i-f-iV) 
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So, finally: 



1 71 

(-i)W ///§VfV ^(-i)" />§Vf s * p ' a±1 " 



4 

With this, (2.19) and 



4 



^-^ re! \ 6 / p/i n , ^-aA r ( 1 



we get 

<i/| e~V |//) 



±(-1)" 



t/H] |//) = lim = 1 (3.20) 



t->o 



(v \e t \u 



if is an integer, else for any \uq) G Ti-kin- Since the result does not depend in 
the choice of \vq), we might take any one. 

From this we see the following: The physical Hilbert space is not separable, it 
consists of the linear forms 

Is := < A I' ^ [0,4/io) 

Ae 5+4/uo Z 
with the physical inner product 

( l& I )phys = ^8e ■ 

So it decomposes into uncountably many one-dimensional, orthogonal subspaces: 

'Hphys = C. 

5e[0,4 M ) 

This shows an important fact: Although the classical system under consideration 
has zero degrees of freedom, the corresponding physical Hilbert space H p h ys of 
the quantum theory has uncountably many dimensions. Also, the theory is not 
superselected, since the operators 

6 X := sin(Ac), A > (3.21) 

are self-adjoint (with respect to the physical inner product) and commute with the 
constraint, hence turn into Dirac observables on the physical Hilbert space. These 
observables act irreducibly in Hphys, hence the physical Hilbert space is not super- 
selected. 



The fact that the quantum theory behaves so significantly different from the 
classical one can be traced back to the quantization procedure. Since the constraint, 
the phase space function c 2 , has no quantization on TLkim one has to come up with 
a different operator and chooses //q 2 sin 2 hqc for some fixed This introduces a 
scale into the theory that has not been there before. Had one chosen the uncountable 
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set of constraints A 2 sin 2 Ac for A > 0, which is classically equivalent to c = 0, then 
the only solution for this system would have been the linear form 

I = J>|- (3.22) 

In this case the space of solutions to the constraints had been one-dimensional right 
from the beginning, but the definition of the physical inner product would have 
been harder (though unnecessary in this simple case). 

Note that the unphysically large size of TLphys should be considered as an artifact 
of the model, resulting from the introduction of the scale hq to a formerly scale- 
invariant theory, i.e. going over from (3.15) to (3.16). This feature does not occur, 
if the model already containes scales, i.e. determined by massive matter fields or the 
cosmological constant A (as, for instance in [11], where the physical Hilbert space 
is onedimensional, as expected). Still, the results, of course, depend on the relation 
of the scales from the theory and the scale from as one can see i.e. in (3.14). 

3.3 The Lorentzian case 

The computation of the physical inner product in the case of vanishing cosmological 
constant has been possible explicitly, because the constraint (3.11) had been changed 
into the classically equivalent (3.12). Unfortunately, this trick is not available for 
the full theory, for some important reason: To become an honest operator in full 
Loop Quantum Gravity, the density weight of a phase space function has to be +1! 
But multiplying with a nonzero power of the triads Ef (which correspond to p in 
the cosmological case) changes the density weight, hence destroys the quantization 
scheme. The honest constraint is therefore 

tf _(*w5*y ^ (3 . 23) 

We would like to consider this constraint instead of (3.16). In this form, the physi- 
cal inner product can no longer be calculated explicitly, but we will show how the 
complexifier coherent states will be able to provide an approximation scheme for it. 

We start by giving all solutions of the constraint equation 

I o C-t = o. (3.24) 
In particular, the solutions are, for 5 € (0,4//o): 

ha,b) = E t; tW, «> 6gC > ( 3 - 25 ) 

whereas the solutions in the sector 5 = are the following: 

ha,b) = E 7; -s-<A| + 6<0|, a,kC. (3.26) 

XeS+A^Z l A + W)P - I A — Mol 2 

Notice that the space of solutions is two-dimensional in each case, consisting of one 
bounded and one unbounded solution. This means that the set of coefficients i( Qj 6) |A) 
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in (3.25) and (3.26) are either bounded or unbounded. The bounded solutions are 
normalizable in the kinematical inner product for the case 5 = only! 

The Master constraint is quite simple, since there is only one constraint: 

M = C^C, (3.27) 

i.e. 



M = /Nsgnp (^^] sgnpv^ (3.28) 
V W / 

From this form it is immediately clear that the solutions (3.25) and (3.26) satisfy 

l(a,b) ° M = 0. (3.29) 
Still, there are more solutions to (3.29). These are, in particular, for 5 € (0, 4/xq): 



kc,d) = E ] + ^ prW, c, deC. (3.30) 

The linear forms (3.30) solve (3.29), but not (3.24). So, as one would have guessed 
from the outset, there are more solutions to the Master Constraint than to the 
original constraint, and it is common belief that the additional solutions will turn 
out to be spurious, that is, have zero norm in the physical inner product yet to 
compute. 

The form of the master constraint (3.28) is too complicated to give a closed 
expression for (u\ M 2n \fi) , which is needed in order to compute the rigging map 
(2.19). It is at this point, where the coherent states come into play, since they allow 
for an approximation of (3.28), as will be shown in the following. 



3.4 Coherent states 

The comp lexifier coherent states have been defined in [7, 8, 9] and their properties 
have been investigated in [17, 20, 21, 18], where their particular use for approxima- 
tions has been pointed out. We will use the simplest of these states, in particular 
the ones for the gauge group U(l). They are defined on each 7is = span (|A) , A G 
5 + 4^o^}, with 5 £ [0, 4// ), by 

\u) = Yj e' x2s ^ e iXu \\). (3.31) 

The complex number u = c — ip that labels the coherent states resemble the phase- 
space coordinates. With our conventions here, the states are peaked in phase-space 
around the points (—2c,p/s 2 ) with width s. Hence, s plays the role of a " classicality 
parameter" . 

The coherent states have a number of properties that make them useful for 
approximations. In particular, they provide a resolution of unity via 
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/T ^rf dpe ^ l c_i P)( c_i Pl = ic W ( 3 - 32 ) 

MO ^ 

Furthermore, the coherent states approximate classical states in the following sense: 
Given a phase space function /(c, p), the expectation value of the operator / := 
f(c, p) in a coherent state |c — ip) is, up to corrections in t, given by 

(c-ipl/lc-iri 
(c — ?p|c — zp) 

where the quality of the approximation depends on /, and will in general not be 
equally good for all phase space functions [18]. 

To calculate the rigging map corresponding to the Master constraint M accord- 
ing to (2.19), we need to calculate 

- 2 00 i / _i \ n 

Me-^M :=£-( — ) (,v\M 2n \v), (3.34) 

which exists, since the basis vectors can be shown to be analytic vectors for M. 
Because of (3.34) we need to know matrix elements of even powers of M in arbitrary 
basis states \i>), This cannot be done analytically, due to the inconvenient 
structure of (3.9) in the operator M, but the coherent states will provide a way of 
approximating (v\M 2n \/j,), in the following way: 

With / = M 2 , M being the master constraint of our model, we will find a 
generalization of (3.33), in particular 

(d - ip'\ f \c - ip) ~( , , (p + p 1 ) + i(c - d) 



(d — ip'\c — ip) 



/(-(■= + C) + <(p -A 61 + ^ (C "^ ). (3-35) 



where / is the analytical continuation of / in both variables. This continuation is 
possible for / = M 2 , since M 2 is in fact real analytic in both variables. Hence, 
to approximate (u\M 2 \/j,), we first compute the matrix elements of M 2 in coherent 
states (u|M 2 |f), with u = c\ — ip±, v = C2 — ip2, and then use (3.32) to write 



HM 2 |/,> = ^ f f d Pl d P2 e- P - 1 ^ (v\v)(v\M 2 \u)(u\ri. 

vrt J\_jl. vli2 2tt 2ir J R 2 

1 MO MO ' 

This will enable us to develop an approximation of (v\M 2 \fi). After that, we will 
employ the completeness of the \u) -basis on 7is to compute 

HM 2 » = HM 2 |A 1 )(A 1 |M 2 |A 2 )---(A n _ 1 |M 2 |/i). (3.36) 

Ai,A 2 ,--- ,A n _ie5+4/^ Z 

Since there is no closed expression for the result available, we will use another ap- 
proximation in order to arrive at a tangible expression for (i^|M 2ri |/i), the quality 
of which will depend on u, fj,, making the approximation phase-space dependent. 
Afterwards, we will plug these approximations for {v\M 2n \n) into (3.34) and inves- 
tigate the limit t — > 0, in order to arrive at a final expression for the rigging map 
(2.19). 
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4 Computation of (u\M z \/j,) 

4.1 Expectation value of the Master constraint in co- 
herent states 

With 

. 3 3 

sgnp VIpI \ v ) = a \ v ) 

= ■ PuW), 

the action of the squared Master constraint on a vector \u) is, via (3.28) given by 
1 



M2 !"' = 256 



PvPv+^oPu+l&fio \v+lQpo) 



-4 PvPu+nm* (pl+8^0 + Pl+4uo) W + 12 Po) 

+ PuPu+8^ {6pl+8v, + 16 PI+4IM> + \ U + 8 P0) 

- Pvpv+An M+8 m + 2A Pt+i w + 2 ^pl + 4 /^-4 W ) W + 4po) 

+ pI(pI+8ho + 16 P^+4 W) + 36 PI + 16 ^-4 W + ^-8 M0 ) H 

M-tipo + 24 pl-4»o + 24 ^ + 4 ^+4 W) ) W~ 4 P0) 

-4 p v p v -n m (pI-8^ + pjL-mo) I" ~ 12 ^°) 

+ Pvpl-B/ioPv-ieiio W - 1Q Po) ■ (4.1) 
The coherent states in the sector 5 £ [0, 4/xo) are given by 

\u) = e- A2 ^e iAu |A). 

Ae5+4^ Z 

As shown in [11], this state, labelled by the complex number u = c — ip, is (in our 
notation) sharply peaked around the point (—2c,p/s 2 ) in phase space with width s. 
One gets, for two coherent states \u) , \v) : 

(v\M 2 \u) = (4.2) 

2 

/ j e 2 e Pii+afi P^i+biJ,oP^+c^( j P^+diJ,o t J u ,fi+2eiM 

terms i/,//e<5+4/ioZ 

with integer numbers a, 6, . . . , e. In order to evaluate this, we write 

Px+atoPx+bwPx+cmtPx+dno = { x + m Vo) 2 + ^o A + #0*0 ( 4 -3) 
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where 

i a + b + c + d 



m 



1 1 

6 ~ 16' 



A = — ^ — -f-(a + b + c + d) 2 - 



a 2 + b 2 + c 2 + d 2 



(4.4) 



and limui^oo xR(x) < oo. If we neglect R(x), we get 

( V \M 2 \u) « ^ e -A 2 S 2 e -(emo) 2 S 2 e iA(«-0) e i(e W) )(«+«) [(A + mfl ) 2 + /UqA] 



terms A 



= ^ e -(^) 2s2 e -(^) 2s2 e -(^— «) 2s2 
terms 

x J] e"( A -^) V [(A + m^o) 2 + ^A] , 
with m = m' — e. Applying the Poisson summation formula we get 



(4.5) 



^ e ( A l ~272~) s (A + m^o) 
Mo 4^ J 



drr ' ( . ^) ° 2 [(x + mfio) 2 + ulA] (4.6) 



1 Hk (■a^.\ 2 s 2 — ^V) 2 s 2 

fl V S 2 



Reinserting this into (4.5) we thus obtain 
(v \M \u) ~ — 4/-^ > e w e v 2s w'v 



.u — v . irn \ 1 



(4.7) 



E 

terms 



g ie/io(«+f)g-(eMo) 2 s 2 



Calculating the sum over the different terms explicitly, taking the prefactors in (4.1) 
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Jefj, (u+v) e -(efi ) 2 s 2 _ e s2 -§^ 



into account, we obtain 

terms 
terms 



z=i(u+v)/M 



sin 2 



D ie^ (M+^') e -(eMo) 2 s 2 m _ q 



terms 



1 ,2^1 / 193 109 

= — e i , , — - — cos 8unz H cos 6unz 

16 z=i(u+«) w \ 48 p 6 P 



487 „ 55 35 

— cos 4/j, z + — cos 2^ - 



16 z=i(u+v)W) v ' 



So 



1 PlF \ -> 2vrn<K . 



( u—v _ 



(4.8) 



e e^i , N sin 8 2 (fi^-i™) + _l_ + 
z=i(u+i;)/L(o \V 2s mo s / 2s^ y 
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z=i(u+v)fi 



C(z) 



4.2 Matrix elements of M 2 with respect to the basis 
vectors 

We use the result from the last section to calculate now arbitrary matrix elements 
of M 2 . For this, we need the resolution of unity provided by the coherent states 
(3.32): 



Mo 



"ITS- 



TV 

/mo ac r 
- 2 ?r Jr 



dp e 7* \c — ip) (c — ip\ = id-^ 



(4.9) 



so we have: 



{u\ M 2 = -4 / I d Pl d P2 e- Pl ^ (v\v)(v\M 2 \u)(u\v) 

TTS Z J\_JL*W ZTT ITT J R 2 



MO ' MO 



(4.10) 



with u = C2 — ip2, v = c\ — ip\. To make life simpler, we will split (4.8) up into 
several parts before inserting and perform the integrations separately. Expanding 
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the first exponential derivative in (4.8) into a power series, we obtain 

(u — v 7rn \ / ^ 

y v ^ s f^o s / 



(v\M z \u) = — W-t > e "o e 



+ 



l AvF s 2m g 2m 



sm z 



Mo V s 2 ^ l m ^ dz 2m \z=(jio(u+v)) 

E-Mi _(ii^_ ™ ) 2 S 2 / 2 X 

1 /IT 2wnU ^^S 2 Lin s 2 



(4.11) 



Mo V s 



16 z=i(w+^)M0 V ' 



We insert now term by term of (4.11) into (4.10). We start with the zeroth order 
term of the derivative expansion: 



Mo 



'7TS U 



dc\ dc2 



2 , 2 
P]+P2 



ZE JL]2 27T 27T 
MO ' M() 1 



x sin 8 ^ (ci + c 2 - i(j> 2 - Pi)) 



(4.12) 



27rm3 ('P2+Pl~ i ( c 2- C l) 1 j Trn \ „2 
V 2^ MO" 2 / 



P2 +Pi - «(c2 - ci) . 7rn 




2s 2 



Mos z 



+ 



2s 2 



Since the exponentials and the sine function are all invariant under the substitution 
ci — ► ci + this integral has the following form: 



MO dc\ S—y 

- n 

Mn n=U 



27T77. \ 2nin6 



MO 



<■"" - / ^/(ci)x(ci), (4.13) 



where 



X(ci) := J] ^^^]( Cl ) e 



ra=0 



(4.14) 



is a sum over characteristic functions. This gives us 



'irs & Jr 2vr J_j<_ 2tt 

mo 



, 2 2 



2 i 2 



M fdcLfndctf dpi dp2 -{*-$) 4, 6 -(m-3) t eWCl e - lMC2 



x sm 



Mo(ci + c 2 - i(p 2 - Pi)) 



(4.15) 



/ P2+Pl-*(c2-ci)\ 2 S 2 

x x(ci) eV 2s2 / 



P2 +Pi - «(C2 - ci)\ 1 



2s 2 



+ 



2s 2 
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Now consider the integration over c\. Notice that the modulus of the integrand is 
dominated by a Gaussian times a polynomial, where the width of the Gaussian is 
proportional to s. So, for small semiclassicality parameter s, the main contribution 
of the integral over c\ will come from the point where the gaussian is peaked, that 
is — c 2 G [— -j 1 , -j-]. At this point, x{~ c 2) = 1- Therefore we have for the integral: 



fi f dci 



r dcx r mo dc 2 f d ^ 

JR 27T 2TT J R 2 



2 e v ' 2 e 



x sin 



Mo(ci + c 2 - z(p 2 - Pi)) 



x e 



P2+Pl~'( c 2~ c l) ^ „2 
, 2? J 8 



P2 + Pi ~ i(02 ~ gi) 

2s 2 



+ 



2s 2 



(4.16) 



X (l + 0(s°°)) 



We omit the (1 + 0(s°°)) in the following calculations. Shifting the integration of 
ci by C2 and rescaling the integration of c 2 by /j,q, we obtain 



'7TS U 



dc\ f n dc2 



6 ./IB 27T 



2vr 



r 7^ 2 R ^ ITj 2 e WC l 1 M0 C 2 ^5- 



l 2 . 2 



dpi (ip2 e 



x sin 



2c 2 +/U (ci - i(p 2 ~Pi)) 



(4.17) 



x e 



/ c 1 -i(p 1 + P2 ) \ 2 

"I 51^ ) s 



1 

2^2 



ci - s(pi +P2) 
2s 2 



Note that the integrand is analytic in ci. We now make a shift of the contour of 
integration of C2 by i(p\ + p 2 ). When doing this, one has to take care about the 
convergence of the integral. As one can check, no problem arises in our case. We 
obtain 



1 r d Cl r dc 2 r _ 

= / / 7T~ / dpidp 2 e 

KS b JR 27T J_ n 2lT J R 2 

(2c 2 + mo(ci + 2zp 2 )) e"3* ( J- - 



x sin 



It is now convenient to carry out the integration over c 2 . We notice that 



/_ 



4 + f 



(4.18) 



(4.19) 



when A" is divisible by four. We get: 



! "~M ! 
(-1) 4 «) — = 



256 



'7TS U 



2tT 



^-dpi dp 2 e 



, 2 .2 



2 . 2 
P1+P2 



x e 2 c i e 4 



& ( 2? ~ h ) ■ (4 ' 20) 
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Integration over c\ gives: 

/DO ITS \ ' 4^ / JR2 

x e- 2 m 2 (^) 2 (4.21) 

We now reinsert the correction that we have omitted since formula (4.17) and are 
left with 



256 



v — fl 

-1)^0 



^ 4 M 



" + /M e^ s2 (1 + 0(0). ( 4 - 22 ) 



For all higher orders of the exponentiated derivative in (4.11) we obtain similar but 
more complicated expressions, that have in common the fact that they stay finite, as 
s tends to 0. The integration over the term containing ( in (4.11) is easily performed 
and yields, up to terms of order 0(s 2 ): 



f4 



dc\ dc2 



16\/ vrs 6 J\ * 71 1 2 27r 27r 

1 Mo ' MO J 



/ 



dpi dp2 e 



., PlVs 2 / P2\ 2 s 2 



pj+pj c 2 a 2 
x e lvCl e~^ C2 e"^" e s a? 



E 2nmi ( P2+Pl- i ( c 2- c l) |_j Tn \ 
e mo ^ mo" 2 / 



2=i(u+u)/io 

2 „ 



. ^0 

32 



193 . 

48 *| 



if—n 

4 M0 



109 . 
1.4 + I" 5 ! 



v — \1 
4 M0 



487 r 
+ "12-^ 



4 M() 



55 . 
,2 + "6 



4 MQ 



,3 



+ 48^ 



4 M0 



(4.23) 



(1 + 0(0 



We now go back to the formula for the resolution of unity in (3.32) and notice that 
the right hand side does not depend on s, so 



id Hs = lim id H , 



(4.24) 



lim 



Mo 



"ITS' 



/mo dc f 
--2- 2tt 7r 

MO 



dp e 7* \c — ip)(c — ip\. 



This allows us to drop all terms of order O(s) in our approximation, and we are left 
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with 



i V — jL 
(-1) 4 «J 

256 V ; 



8 

4 + - A 



v + n 



32 V ; 



4 M() 

193 

"48" tel' 4 



+ 



+ 



487 
~12 



U — fJ, 

4 M0 



55 . 



4 M0 



109 

O I 4 M0 I > d 

35 

1,1 + 48 1^1,0 



4 M0 



(4.25) 



Note that the approximation now comes from the neglection of the R(x) - term in 
formula (4.3), but not any more from s-corrections, since we have gone over to the 
limit s — > 0, in which all these terms vanish. 



5 Computation of (is\M 2n \fi) for n > 1 and rig- 
ging map 

5.1 Expectation values of M 2n 

The approximate formula for the matrix elements of M 2 have been calculated in 
the last section. As demanded the approximation vanishes if is not an integer. 
If it is, the matrix elements can be approximated by 



{v\M 2 \fi) « 

(-1)" 
256 



\ 2 

v + H 



A + N 

193 2 436 2x 974 2x 220 2 35 2x 

H ^-Vo°\n\,4 H ^-^o<V|,3 H 5-/^0^1^1,2 H ^-/VV|,i + -t-Mo"|JV|,o 



with JV = ig. So 

("I* 2 ^ = 1 W(4 + 8 W )(^) 2 + ^"-")' («) 
with i? being correction terms that we neglect for the time being. 

Since the \v) with v G S + 4/j,qZ provide a basis of Hs, we can easily compute the 
matrix elements of M 2n : 

(—\\ N / 8 \ / 8 \ f 8 

( ^ 2n|rf " W £ ( 4+^1 ) ( 4 + ^ J"-- 4+^^ 



Ai G <5 + 



X 



4/io / \ 4/i / \ 4/io 



^ + Ai\ 2 /Ai + A 2 \ 2 /A n _i + /U x2 
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It is very difficult to give a closed expression for this formula. Still, we want to show 
a way of calculating these matrix elements in principle. For this, we first rewrite 
the summation: 



(u\M 2n \fi) 



(-1) 4 «> 
256™ 



e n 

Ii,...,A„6Z fc=l L 
Efc 'fe = 



4 + 4 



(5.2) 



n 



V + jU 



fc-1 



j'=i 



i=fc+i 



Note that, whenever is not an integer, the sum is empty and thus, by convention, 
0, as demanded. Because of 



V + jJL 



k— 1 n 

^0^2 1 j + no l i 

3=1 j=k+l 



da\ 



a k =0 



(5.3) 



and 



r 27T db 
Jo 2^' 



(5.4) 



We can rewrite (5.2) in the following way, after shifting each of the sums over 1% by 
-4: 



256 r ' 



2tt 



ii,.,i»a'' u 



^/>fe"E fe ^+4r 

2vr 



n 

,fc=i 



4 



n 

fc=i 



dai I „ 

K afe=0 



a k {^ - MO E.ti 1 »i + MO E?=fc+i «j + 4 M o(n-2fc+l)) 



(-1) ^0 



E 



2tt 



256™ ^ 7 2vr 

l\,...,ln& 



n 

fe=i 



(5.5) 



n <i 2 

n ^2 

fc=l ^ 



a k =0 



n 

fc=i 



da 2 . 



= 4o fe /i (n-2fc+l) JkW) (E*=i a i ~ E?=fc+i a j) 



a k =0 



By the binomial formula, this is 



256 r ' 



2vr 



d 2 

n jo 2 



a k =0 



I] e 4 ^«0(n-2 fc+ l) (1 + a(ai) ^ b)) 8 (5 6) 



fc=l 
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with 

C k ( ai ,...,a n ,b) = ^(E^-E^+i (57) 
It is now straightforward, though tedious, do perform the 2n differentiations. Since 



n 



sgn(l - k) ((1 + Ci(a, b)) n - (1 + Ci(a, b)) 11 - 1 ) , (5 



it is predictable that the differentiations and setting of a = afterwards will yield 
the result 



n d 2 

fc=i aa * 



a k =0 



I] e 4a ^°("- 2fe+1 ) (1 + C fe (ai,...,a n ,6)) 6 



(5.9) 



fc=i 



8n 



^P nim (i/, M )(l+e- i6 ) 8n_m e^, 



m=0 



where the P ntm (v, /x) are polynomials in f and |tx. Since 



/" 27r * e Ki^+ 4 ")(i +e -^-"v<™ = ( 



8n — m 



4n+ 4mo 



m 



and 



8n — m 



4^o 



8n 



8n- 1 



(5.10) 



(5.11) 



' 4n + ^-m + r 
8n — m + 1 



8n 
4 - + g 



the form of the result is known: 



( 1 4 M() / Or, 



(5.12) 



where Q n {y,n) is another polynomial in i/ and /i, consisting of linear combinations 
of the P njm (u,fi). However, it is a considerable computational task to obtain the 
polynomials P n: m(v, fi) for higher m, and, to the best of the authors' knowledge, 
no closed formula for them exists. So, although one can compute the approximate 
matrix elements of (v\ M 2n explicitly order by order, no closed formula for them 
is known. 



5.2 Further approximations 

In order to compute the rigging map approximatively, we start again from formula 
(5.2) 



u — (j, 

-1) 4«T 
256™ 



e n 

h,...,i„ez fc=i 

n 

fe=i 



4 + Z fc 



v + \i 



(5.13) 



A fc (Z) 
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with 



k— 1 n 

Ajfc(f) := MoX^J - E l r ( 5 - 14 ) 

j=i j=fc+i 

Now assume that v — \i is small compared to ^ + /x. Then the summation of the 
Ik ranges over —4, . . . , 4, and since the binomial prefactors are nearly Gaussian 
distributions in the Ik, it is safe to assume that the largest contribution to the sum 
will be made by products, where all the Ik are close to 0, compared to 4, hence the 
Afc(Z) are much smaller that ^j^. So, as a test, we could expand formula (5.13) 
into powers of the sums over the Afc(Z): 



{u\M In \n) 



(-1) 4 "o 
256™ 



v + fi 



e n 

h,...,i„ez k=i L 
Ek 'fe = i^f 



2ra n 

- 2^A fc (f) 
fc=i 

\ 2 



4 + / fc 



+ 2 K> fe (f) - ^Afc(f) 2 



vfc=l 



fc=l 



2n-l 



f + n 



(5.15) 



2n-2 



+ 



For what follows, we assume n > 1, since only then are the Afc(Z) / 0. The cases 
n = and n = 1 are trivial and already known respectively. 
Now first we observe that 



E /« = E i E /(*('"))■ 



n! 

(l,...,in6Z o-eSym„ 

U — fJ, 

4^0 



(5.16) 



where Sym n is the permutation group in n elements and cr(Z) means the vector, 
of which the n entries in the vector I are permuted according to a. Our second 
observation is that 



E A *Hft = °. 

(reSym n fc=l 



(5.17) 



so the "linear" term in (5.15) vanishes, and, keeping only the first two terms of our 
expansion, we are left with 



(v\M 2n \v) 



(-1) 
256™ 



e n 

h,...,ln 6Z fc=1 L 
4mo 



4 + / fc 



Efe 'fc — 4u n 



f + A* 



2n 



+ 2 E A ^( f ) - E A ^ f ) 5 



s.fc=l 



fe=l 



(5.18) 



2n-2 
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(-1) 



256 r 



8n 
4n+^ 



v + n 



2n 



+ 



■1) 4 w / v + n 



256™ra! 



2ra-2 



e n 



h,...,i„ez 

Efc 'fe = 4^ 



fc=l 



4 + 4 



E x>*(*fl) ( 5 - 19 ) 



(TgSym n \fc=l 



e n 
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Let us consider the middle term in formula (5.19). First we note that 



1 " f (ra-l)(/ n -/i) + (n-3)(/ n _i-Z 2 ) + --- + 2(/n + 3 - U-i) n odd 
-VA t (I) = <^ v — — 7 

Mo tZl" ! (n -!)(/„- h) + (n - 3)(Z„_i - Z 2 ) H h(Zn + i-Z«) n even 
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(5.20) 



So, squaring this, one obtains a sum over two different types of terms, in particular 
(l a -l b ) 2 and {l a -l b ){l c -ld), (5-21) 

where in the first term a and b are different and in the second term all a, b, c d are 
different. Just because of this, terms of the last kind cancel when summing over 
Sym„: 



E ( l <r(a) ~ l<r(b))(l<T(c) ~ l a(d)) = 0. 
o-GSym n 

So we can only keep the terms of the first kind and, with (5.20), we obtain: 
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E ((n " l) 2 (/a(n) " W)) 2 + ( n " 3 ) 2 (^(n-l) " ^(2)) 2 + ■ • •) 
creSym n 



= 2A W x 2 1 (n-2)! ^^(/ fc -/ m ) 2 , 



(5.23) 
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l 2 + 3 2 + ■ ■ ■ + (n - l) 2 n even 
2 2 + 4 2 + • • • + (n - l) 2 n odd 



(5.24) 
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A short analysis reveals that k n = |(n — l)n(n + 1), independently of whether n is 
even or odd. So, combining the results from (5.19) and (5.23), we get 
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Let us now evaluate the last line in this expression. With a bit of algebra one can 
show (for n > 1): 

* e E -(»-') (E g) + | ( " + T ' 1)2 f E w,) ■ 

a& ym n fc=l \fc=l / ^ b \ fc <J / 



(5.28) 



of which we only keep the first term, because the second one is small compared to 
it for large n. 
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The constrained sum can be evaluated by means of the same trick that has been 
applied to calculate the constrained sum (5.2). Without showing all the steps again, 
one obtains: 
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With formula (5.11) one can calculate 
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Inserting (5.30) and (5.31) into (5.27) gives: 
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In our calculations we have assumed n > 1, but formula (5.32) also produces the 
right result for n = 1, as one can readily see. Note that the binomial coefficient 
for n = produces a for v ^ [i. For v = \x and n = the binomial coefficient 
gives a 1, as does the whole expression (5.32). So for n = the expression (5.32) 
equals 5 Ufl , as expected. Hence, the approximation formula is valid even for all n > 0. 



5.3 The physical inner product 

We see that, considering v — \t being small compared to v + /x, the first term in 
(5.32) is large compared to the second and the third, although not uniformly for 
all n. Since, to compute the rigging map, one has to sum over all the orders of n, 
the contribution of the two correction terms could be considerable, if v + \i is too 
small. So we see that for small v + [/,, the major contribution might not come from 
the zeroth order term, but from the corrections. Fortunately, one can calculate the 
corrections to the rigging map order by order, which we will do now. The zeroth 
order term is 



The rigging map is given by : 



f — fi 
-1)4^ 

256™ 
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4 - + ^ 
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(5.33) 
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(5.34) 



where \vq) is a reference vector. With the abbreviations 



v + n 



32/io\/i 

we can make use of the integral (3.19) 
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where the integration contour is chosen to pass slightly left from the pole at 0. We 
now apply a coordinate transformation 8s + 1 = — r, so we get 



^ t-ioo—e 



ds r(-s)a 2s 
1 



r(8s + i) 



I6iri 



r(4s + i + iv)r(4s + i-N) 

ioo—l+e 



(5.37) 



r+l 

drT(—r)a 4 



ioo — l+e 



r^-I + iv)r(i-|-iv)- 



29 



Observing that the integrand has no poles in [— 1 + e, —e] x iR, we can shift the 
integration curve and obtain 



16vri J_ ic 



— r + 1 

dr r(— r) a 4 



r(^-| + ^)r(i-|-iv) 



« a 4 2^ 
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(5.39) 



This means that, reinserting a and N from (5.35), only the n = O-Term survives in 
the limit: 
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Together with 



(-l) 4 ^o 



(5.40) 



and (5.34) we arrive at the result 
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The above calculation can be carried through likewise with the correction terms in 
(5.32), which yields, after some careful analysis: 
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For the calculations with the corrections, the choice of uq results in a choice of 
normalization of the rigging map. However, the dependence of the normalization 
on uq is complicated. In (5.43) a particular normalization has been chosen, to make 
the formula comparable to (5.42). 



5.4 Discussion 

First note that the choice of the reference vector \vq) is only a choice of normalization 
of our rigging map. It does not even have to be from the same sector 5 6 [0,4/xq) 
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than v and \i. Note further that in our approximation (5.15) we assumed that v + fx 
is large compared to v — fi. In particular, the zeroth order approximation (5.42) 
cannot be trusted if this condition is violated and, as one can see, it even diverges if 
/j, approaches —v. We will now discuss this formula for the sectors 5 = and (5/0 
separately. 

5 € (0, 4/i ): 

The rigging map we have produced suggests, as a zeroth order approximation, 
that 

« « E <A|, (5.44) 

where the result should only be trusted, where A « v, and both A and v are large. 
We know that by construction (2.19), solves the master constraint, i.e. 

r][\v)] o M = 0. (5.45) 

In chapter 3.1 we calculated the solutions to the constraint C explicitly. There we 
saw that in each sector 5 £ (0, 4/io) there are two solutions (3.25), one of which 

is bounded and one of which is unbounded. For all sectors 6^0, the unbounded 

i i 

solution goes like |A|5, whereas the bounded one goes like | A | a for |A| — > oo. 
Furthermore, the master constraint M has two further solutions in each sector, 
given by (3.30), which behave asymptotically like |A|a and |A|5, respectively. 

For fixed \ the approximation stays bounded in the regime where it is supposed 
to be a good approximation to the real and goes asymptotically like |A| _ 4. 

Although this is not quite the right asymptotical behaviour, it is safe to assume 
that the true r}[\v)] is the bounded solution 

« = E ^ pr^-r ( 5 - 46 ) 

i.e. the one going like | A| 2 . So, although the approximation (5.44) does not recover 
the exact asymptotical behaviour of the physical solution (5.46), it is good enough 
to identify the right solution (5.46) among the four possible ones (3.25) and (3.30) 
in each sector S 6 (0, 4/xo). 

Numerical comparison of the different approximations (5.2) and (5.32) with the 
real values for (u\M 2n \fi) for different n indicates that the approximation provided 
solely by the coherent states (5.2) is quite good, whereas (5.32) is not quite that 
accurate. This shows that the reason for the slightly different asymptotic behaviour 
of the final approximation of the physical solution compared to the one expected is 
caused by further simplifying the formula in chapter 5.2, in particular by by going 
over from (5.13) to (5.18). This was necessary in order to bring the expressions 
into a form such that further calculations were possible, but it happened at the 
cost of accuracy. Still, the approximation formula can in principle be improved by 
computing and adding further corrections, such as happened with the first order 
corrections in (5.43). 
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Note that in particular, although the rigging map was defined by the direct inte- 
gral decomposition of M (2.19), the physical Hilbert space contains only solutions 
of M that are also solutions of C. This shows that, although one had replaced C 
by M, no information is lost, and hence the master constraint programme proved 
to be the right tool to define the physical Hilbert space in the presence of a non- 
self- adjoint constraint. 

Note further that since, analogously to the case of chapter 3.2, in each sector 
5 £ (0, 4/xo) there is exactly one solution of nonzero physical norm (hence all other 
solutions in the same sector are spurious), the Hilbert space is again non-separable, 
in particular it contains 

C C H phys . (5.47) 

56(0,4^0) 

Hence, in the sectors 5 G (0, 4/xo), the result for the physical Hilbert space provided 
by the master constraint is the same as for the simplified self- adjoint constraint 
(3.16), which in fact is classically equivalent to (3.23) apart from p = 0. So the 
more complicated, non-self-adjoint constraint (3.23), together with the master con- 
straint leads to the same physical Hilbert space as the simplified constraint (3.16). 

5 = 0: 

The case of 5 = has to be treated differently than the case 5 £ (0, 4^o)- The 
reason for this is that the vector |0) is contained in Ti.s=o, which corresponds by 
(3.6) to the classical singularity p = 0, defining the points in phase space where the 
two constraints (3.16) and (3.23) differ classically. Indeed, |0) is an eigenvector with 
eigenvalue of (3.23) 

C\0) = ( s -^.) 2 sgnp^pt|0) = 0, (5.48) 
\ W / 

but not of (3.16). Hence the solutions of (3.23), given by (3.26), have a different 
structure than (3.25): There is no bounded solution with asymptotical behaviour 
|A|~2 5 instead there is the eigenvector |0). This eigenvector constitutes TCq P and 
has to be excluded from the following discussion, which deals with 7{q c only (see 
chapter 2.3). 

Although there is no analogue of (5.46) in the sector with 5 = 0, the approxi- 
mation that has been calculated gives 

viw)} « E -^—r <A|, 

even for v € 4/ioZ, i.e. in the sector 5 = 0. So the approximation scheme treats 
all 5 £ [0, 4//o) equally. In the case for 5^0, (5.44) has been interpreted as an 
approximation to the solution (5.46), which, as we have seen, does not exist in this 
sector. Thus, it is not clear which solution the rigging map approximates, if any. 
We attempt to interpret this as follows: 

The disappearance of an analogue to the solution (5.46) for 5 = can be traced 

back to the appearance of sgnp-y/jp[ in the constraint (3.23), which results in an 
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eigenvector of C (hence M) with eigenvalue 0. The eigenvector |0) of C has to 
be removed from the analysis, since discrete and continuous spectrum have to be 
treated separately in the direct integral decomposition 2 . Hence, the Hilbert space 
Hq has to be replaced by span{|z/} | v 6 4 / uo^\{0}}, which is a sub-Hilbert space 
in Hq of codimension 1. The coherent states, however, are defined on all of Ho, 
in particular the projection to |0) is nonzero for all coherent states in Ho, i-e. 
states (3.31) with 5 = 0. So, by using the coherent states for the semiclassical 
approximation on the sector 5 = 0, one is extending the space one is allowed to 
look for solutions. So we have to conclude that the approximation provided by the 
coherent states are not appropriate for the sector 5 = 0, since, strictly speaking, 
there are no coherent states living on this sector. Thus, the approximation cannot 
be trusted at all, and we conclude that there is no physical solution in the sector 

5 = 0. Hence, the whole physical Hilbert space is 

H phys = Hl P @H a o c = span{|0)} © C. = C (5.49) 

56(0, 4 M0 ) <5e[0,4 Mo ) 

6 Summary and conclusion 

In this article, we computed the physical inner product of a simple cosmological 
model. We did so in two different ways: Firstly, we rewrote the classical side before 
quantization, in order to arrive at a much simpler constraint on the quantum side, 
the physical inner product of which could be computed easily. In each Hilbert space 
Hs, there are a bounded and an unbounded solution to the constraint, and the 
bounded one was constituting the physical Hilbert space H p h ys , while the other 
solution was "spurious", i.e. not contained in it. In the model calculated H p h ys is 
non-separable. This unphysical feature is a result of the quantization procedure, 
which brakes scale-invariance of the classical theory. Still, this was no particular 
problem for our calculations. 

Secondly, we attempted to compute the physical inner product directly, with- 
out rewriting of the classical side. This resulted in a system which is classically 
equivalent up to points with p = on the phase space. The constraint of the new 
system is non-symmetric, which is why we employed the master constraint method 
to arrive at a self-adjoint constraint. Since the formulae for this operator were too 
complicated to compute exactly, we relied on the complexifier coherent states to 
approximate the physical inner product. This has led to interesting results: 

Similarly to the case already observed, on each sector 5 G (0, 4/io) there are a 
bounded and an unbounded solution, showing that rewriting on the classical side did 
not change the quantum theory significantly on the parts of the phase space where 
the two classical theories are equivalent. The sector 5 = contained the eigenvector 
of p belonging to vanishing momentum p = 0, which defines the phase-space points 
where the two classical theories do not agree, and indeed the solution structure of 
the quantum constraints changed significantly. What did not change is the physical 
Hilbert space. Although now containing an eigenvector, |0), the resulting physical 
Hilbert space Hphy S was isomorphic to the one computed from the simpler model. 

2 In particular, one has to perform the direct integral decomposition separately for each part of the 
spectrum, leading to a different choice of normalization vector \v ) in (2.19) 



33 



So rewriting the classical theory changed the structure of the solution space in 
quantum theory at the particular points, where the classical theories disagree, but, 
after going over to the physical Hilbert space, leading to the right result. 

The approximation of the physical inner product with the help of the complexifier 
coherent states was good enough to indicate that for each sector 5 G (0,4^ ) the 
bounded solution was chosen to constitute the physical inner product, while all 
the unbounded solutions turned out to be spurious. This demonstrated firstly that 
going over to the master constraint did not change the quantum behaviour. In 
particular, although there are many more solutions to the master constraint M 
than to the original constraint C, computing the physical inner product with M 
still resulted in a physical Hilbert space consisting of solutions to C only. Hence 
the additional solutions provided by the master constraint did not contribute to the 
physical Hilbert space, showing the strength of the master constraint programme. 

This demonstrates that the complexifier coherent states provide a way to approx- 
imate physical inner products. Furthermore, the corrections to the approximated 
rigging map can in principle be calculated order by order (the "perturbation pa- 
rameter" in our case being the distance to the region of phase space where the 
approximation was to hold best) and be added to the previously obtained result, 
granting a way to improve the approximation by adding more and more correction 
terms. 
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